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I.  INTRODUCTION. 

Hydrogen  la  the  most  important  high-energy  fuel  in 
rocket  technology.  Hydrogen-oxygen  rocket  motors  are 
already  being  utilized,  and  other  liquid  propellant 
systems  using  hydrogen  with  several  oxidizers  are  in  a 
well  advanced  state  of  development. 

Hydrogen  has  very  distinctive  physico-chemical 
characteristics.  It  has  unusually  low  critical  pressure 
of  the  order  of  12  atmospheres  with  the  result  that  it  is 
injected  into  the  c  >mbustion  chambers  of  rocket  orators  at 
supercritical  conditions.  At  pressures  higher  than  the 
critical  value  the  surface  tension  becomes  equal  to  zero, 
and  therefore  hydrogen  cannot  form  droplets  within  the 
combustion  chamber,  acting  as  a  homogeneous  fluid. 

On  the  other  hand,  most  oxidizers  have  critical 
pressures  well  above  the  normal  working  pressures  of 
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liquid  propellant  rocket  motors.  Therefore,  the  basic 
combustion  process  in  a  rocket  motor  occurs  as  burning 
of  a  spray  of  oxidizer  droplets  within  a  hydrogen  atmos¬ 
phere. 

The  object  of  this  paper  is  the  study  of  the  basic 
process  of  combustion  of  single  oxidizer  droplets  in  a 
quiescent  hydrogen  atmosphere. 

There  exists  an  abundant  literature  on  droplet 
combustion  (refs.  1,2,3  and  4,  for  example).  Most 
studies  refer  to  the  case  of  fuel  droplets  burning  in 
air.  Some  studies  consider  more  sophisticated  fuels 
and  there  a re  a  few  studies  on  monopropellant  droplets. 

On  the  other  hand,  no  specific  theoretical  or 
experimental  studies  have  been  performed  on  the  combustion 
of  oxidizer  droplets  in  hydrogen.  This  case  cannot  be 
included  in  a  general  theoretical  study  on  droplet  combus¬ 
tion  owing  to  the  specific  properties  of  the  hydrogen. 

Its  low  density  and  low  molecular  weight  make  diffusion 
conditions  especially  critical  and  the  assumption  of 
taking  average  values  for  gas  density  and  for  transport 
eoeficients  regardless  of  mixture  composlticn,  which  is 
normally  adapted  in  droplet  combustion  studies,  may  in  - 
troduoe  important  errors  for  the  case  of  hydrogen  and  it 
is  shown  in  the  paper. 


Some  results  on  the  combustion  of  nitric  acid  and 


bromine  droplets  in  hydrogen  were  published  by  ourselves 
in  Refs.  4  and  5»  These  studies  have  been  completed  and 
combustion  of  oxygen  droplets  has  been  also  included  • 

Bromine  was  selected  because  chemical  kinetics  of 
the  H^-Br^  reaction  is  well  known,  and  this  offered  a 
good  opportunity  for  the  study  of  the  influence  of 
chemical  kinetics  on  droplet  combustion. 

Nitric  acid  was  chosen  because  this  liquid  oxidizer 
is  particularly  suitable  to  carry  out  experimental  studies 
and,  finally,  oxygen  was  selected  because  of  its  high 
technical  interest.  This  last  study  is  an  application  of 
a  general  study  on  the  hydrogen-oxygen  diffusion  flame  with 
spherical  symmetry,  shown  in  Refs.  6  and  7* 

2  “  MODEL  OF  THE  PROCESS 

Combustion  of  oxidizer  droplets  in  hydrogen  atmosphere 
will  be  studied  under  the  following  assumptions  : 

1  -  Spherical  Symmetry.  Therefore,  natural  convection 

effects,  always  small  in  a  hydrogen  atmosphere, 
will  be  disregarded. 

2  -  Stationary  conditions.  Recession  velocity  of  the 

droplet  surface  will  be  disregarded  as  compared  with 
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diffusion  velocities  oi'  chemical  species,  and  the 
energy  and  time  required  to  heat  up  the  droplet  up 
to  its  boiling  temperature  will  be  also  disregarded. 

3  -  Constant  pressure  since  it  is  practically  constant 

in  this  type  of  processes. 

4  -  Three  chemical  species:  Hydrogen,  oxidizer  and  reac¬ 

tion  products  will  only  be  considered.  The  influen¬ 
ce  of  radicals  will  only  appear  through  the  chemical 
kinetics  of  the  process. 

5  -  An  over-all  reaction  rate  will  be  taken  in  order  to 

consider  chemical  kinetics  of  the  process.  This 
over-all  reaction  rate  may  be  derived  from  the  complete 
chemical  reaction  model  of  the  combustion  process. 

6  -  A  linear  variation  of  gas  density  and  thermal  conduc¬ 

tivity  as  function  of  mixture  composition  will  be 
considered. 

3  -  FUNDAMENTAL  EQUATIONS 

The  fundamental  equations  of  the  process  will  be 
expressed  in  terms  of  the  dimensionless  variables  : 

c°  1 

4  n  1 


x 


r 


(1) 
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c° 

0  =  (T  -  T  )  (2) 

o  o 

and  the  over-all  combustion  reaction  will  be  assumed  to 
be  of  the  form  : 


M2  =  $2  M3 


(3) 


Under  the  aforementioned  assumption  and  in  terms  of 
x  and  0  the  fundamental  equations  of  the  process  are  as 
follows : 


a)  Continuity 


d& 


i 


dx 


(4) 


which  expresses  the  conservation  of  mass  of  each 
chemical  species  i  • 


b)  Energy 


(  6j“l  )  +0^1  +  6^-1 ) 


,  A  d0  +  0  o 

f  dx 


(5) 


This  equation  expresses  that  the  sum  of  the  energy 
fluxes  of  heat  and  enthalpy  is  constant.  This  constant  has 
been  obtained  by  expressing  at  the  droplet's  surface  that 
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the  heat  received  by  the  droplet  is  utilized  in  evaporating  the 
liquid. 

In  this  equation  parameter  is  : 
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heat  of  reaction  n° 

*r 

is  given 

by 
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c )  Diffusion 
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Assuming  that  only  concentration  diffusion  exists,  it 
results  : 
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in  which,  oCiJ  is  the  Lewis -Semenov  number  given  by  : 
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ij  • 


,.«i  _x ... 
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(9) 


d)  Boundary  Conditions 

Boundary  conditions  arc  as  follow*:  : 


At  th«  droplet's  surface  : 


x 


0 


At  the  infinity  : 


x 


(10) 


(11) 


4  -  SOLUTION  OF  THE  SYSTEM 

In  the  first  place  the  problem  will  be  solved  by  assuming 
an  infinitely  fast  reaction  rate,  which  is  equivalent  of 
taking  a  zero- thickness  reaction  zone. 

This  is  the  well  known  Burke  and  Schumann  assumption 

4 

for  studying  flames  and  it  has  been  proven  that  it  gives  a 
good  approxlsiation  of  results  in  droplet  combustion  provided 
that  the  droplets  are  not  too  ssiall  or  the  pressure  is  not 
too  low. 

In  this  case  in  which  hydrogen  is  the  reactant  atmosphere 
with  a  molecular  weight  much  smaller  than  that  of  the  oxidizer 
droplet,  it  was  considered  that  the  influence  of  mixture 
composition  on  the  value  of  density  and  thermal  conductivity 
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might  be  specially  important.  Therefore,  linaar  variationa 
of  danaity  and  tharmal  conductivity  aa  functions  of  mass 
fraction  of  tha  combustion  products  wars  takan: 


_JL 

Txi 


_A_ 

I 


ci3  *  gi3  Y3 


(12) 

(13) 


With  tha  Burks -Schumann  assumption,  tha  problasi  has 
an  axact  analytical  solution. 

Mass  fraction  and  tampers tura  ara  given  by  : 


ia  obtained  by  expreaaing  in  tha  aforegoing  aquationa  that 
for  x  b  ia  O  *  0(  . 

Flam*  front  poaition  x  =  x*  and  maximum  temperature 
0*  are  obtained  by  expreaaing  that  for  x  =  x*  ia  =  1 

From  the  value  of  xfl  the  evaporation  conatant  k, 

2  #  x 

equal  to  the  alope  of  curve  rg  =  f(t)  ia  obtained  by 


f 
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means  of  the  expression  : 

dr^  A 

k . 5—  =  -  -----  (19) 

dt  4 

5  -  THEORETICAL  RESULTS  FOR  INFINITE  REACTION  RATE 

Theoretical  result*  of  the  dimensionless  burning  rate 
f  evaporation  constant  k  and  flame/droplet  radius  ratio 

*.•/ r 

a  were  obtained  for  nitric  acid,  bromine  and  oxygen 
droplets  assuming  an  infinitely  fast  reaction  rate. 

Results  are  shown  in  Table  No.  1,  in  which,  for  com¬ 
parison,  results  corresponding  to  the  combustion  of  hydro¬ 
carbons  droplets  in  air  are  also  included,  as  well  as  some 
other  representative  propellant  combinations. 

It  ma.v  be  seen  that  combustion  in  hydrogen  is  con  - 
siderably  faster  and  that  the  flame  is  much  closer  to  the 
droplet  surface  than  in  the  other  propellant  combinations  . 

It  is  also  shown  in  the  Table  the  very  important  in¬ 
fluence,  for  the  case  of  combustion  of  oxygen  droplets  in 
hydrogen,  of  considering  constant  or  variables  values  of 
the  gas  density  and  transport  coefficients  as  function  of 
mixture  composition. 
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6  -  FINITE  REACTION  RATE 

System  of  differential  equations  (4),  (5)  and  (6)  was 
solved  for  the  case  of  finite  reaction  rate  by  utilizing  the 
approximated  analytical  method  described  in  Refs.  4  and  5  • 

For  applying  this  analytical  method  an  over-all  reaction 
rate  is  needed.  Therefore,  it  could  not  be  utilized  for  the 
cases  of  nitric  acid  and  oxygen  for  which  such  over-all 
reaction  rates  were  not  available. 

For  the  bromine-hydrogen  reaction  the  Campbell  and 
Hirschfelder  reaction  model  was  adopted  from  which  an  over¬ 
all  reaction  rate  was  obtained  (see  Ref.  5)* 

Results  were  obtained  as  function  of  the  dimensionless 
product  of  the  3/4  th  power  of  the  pressure  multiplied  by 
the  droplet  radius  in  the  form  : 


x  • 

s 


¥><  Ql/2  P5/4 


(20) 


x 


y« l1/z  p3/4  r  ) 
8 


(21) 


in  which  parameter 


0  is  given  by: 
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Q  =  1.8816  109  M,  (M_/M,)1/2  (UT_)"  3/2  C  3  _1  (22) 

j  2  1  1  p 

The  values  of  x  and  r*  /r  are  shown  in  Figs.l 
and  2  in  which  the  values  for  nitric  acid  and  oxygen  dro¬ 
plets  for  the  case  of  an  infinitely  fast  reaction  rate 
are  also  shown. 

Results  show  that  for  the  case  of  bromine  droplets 
combustion  does  not  exist  under  a  minimum  value  of  either 
the  pressure  or  the  droplet  radius*  .  it  may  also  be  seen 
that  above  this  value  results  tend  rapidly  towards  those 
obtained  assuming  an  infinitely  fast  reaction  rate. 

7  -  EXPERIMENTAL  RESULTS 

Experimental  results  were  obtained  by  burning  droplets 
suspended  of  a  thin  quartz  fiber  in  a  chamber  filled  with 
hydrogen  at  different  pressures.  Photographic  records  were 
taken  of  the  droplets  diameter  and  flame  diameter  as  func¬ 
tion  of  time 

•  .-  This  conclusion  was  established  in  general  for 

droplets  combustion  considering  finite  chemical  kinetics 
in  Ref.  4  • 


Dimensionless  burning  rote 


FIG.  1 
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These  •xporimental  studies  vert  only  carried  out 
for  nitric  acid  and  bromine  droplets  since  it  was  not 
possible  to  carry  out  this  experimental  work  for  oxygen 
droplets. 

Experimental  results  for  nitric  acid  droplets  are 
shown  in  Fig.  3«  and  numerical  results  obtained  in  Table 
No.  1,  in  which  it  may  be  seen  that  there  is  a  reasonable 
agreement  between  theoretical  and  experimental  results. 

Combustion  of  bromine  droplets  in  hydrogen  is  of 
an  explosive  type,  as  high-speed  photographs  show  in  Fig. 

4.  A  similar  type  of  combustion  has  been  observed 
sometimes  in  metal  particles  and  there  not  yet  exists  a  clear 
understanding  of  this  phenomenon. 

Due  to  this  type  of  rough  combustion  comparison 
between  theoretical  and  experimental  results  was  not, 
lor  this  case,  carried  out. 

8  -  CONCLUSIONS 

a)  Coisustion  of  liquid  oxidizers  in  hydrogen 
ocours  as  burning  of  a  spray  of  oxidizers  droplets  in 
an  homogeneous  hydrogen  atmosphere. 


square 
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•5%  Nitric  acid  -» - f- k  =  5.8KlO',cmV*«c 

98%  Nitric  acid  ft  *  6.5 *10**  cm*/sac 

98%  N.a.(Th«or«tical) - R«4.4  *10'*  cmVstc 
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BURNING  HATES  ,  EVAPORATION  CONSTANTS  AND  KLAME-UKOP1-.1ST 
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Data  taken  from  our  own  experimental  results 
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b)  Combustion  of  oxidizer  droplets  in  hydrogen  is 
much  faster  than  combustion  of  most  usual  combinations  of 
propellants . 

c)  The  hydrogen-oxidizor  flame  is  located  very 
close  to  the  oxidizer  droplet  surface. 

d)  When  calculating  combustion  of  liquid  oxidizers 
in  hydrogen  it  is  important  to  take  into  account  proper 
values  of  gas  density  and  transport  coefficients  as  func¬ 
tion  of  mixture  composition. 


LIST  OF  PRINCIPAL  SYMBOLS 


C®  -  molecular  heat  capacity  at  constant  prasaura  (constant)* 
D±j  -  diffusion  coafficlant 
A  -  mass  flow* 

M  -  molecular  weight 

p  -  pressure 

q°  -  mole]  heat  of  evaporation  of  liquid  oxidizer 

q°  -  molal  heat  of  reaction 

r  -  radius 

w^  -  reaction  rate* 

x  -  dimensionless  coordinate 

-  mass  fraction  of  species  i 
6^  -  flux  fraction  of  species  i 

■J,  -  average  value  of  thermal  conductivity 

*)  -  stoichiometric  coefficients* 

Tij  ”  average  value  of  density  of  species  i  and  j  . 

0  -  dimensionless  tesq>erature 

Subscripts 

1,2,3  -  oxidissr,  hydrogen  and  reaction  products,  respec¬ 
tively* 

s  -  at  the  droplet  surface. 
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The  object  of  this  paper  is  the  study  of  the  basic  process  of  combustion  of  single 
oxidizer  droplets  in  a  quiescent  hydrogen  atmosphere.  There  exists  an  abundant 
literature  on  droplet  combustion.  Most  studies  refer  to  the  ctse  of  fuel  droplets 
burning  in  air.  Some  studies  consider  more  sophisticated  fuels  and  there  are  a  few 
studies  on  monopropellant  droplets.  On  the  other  hand,  no  specific  theoretical  or 
experimental  studies  have  been  performed  on  the  combustion  of  oxidizer  droplets  in 
hydrogen.  This  case  cannot  be  included  in  a  general  theoretical  study  on  droplet 
combustion  owing  to  the  specific  properties  of  the  hydrogen.  Its  low  density  and 
low  molecular  weight  make  diffusion  conditions  especially  critical  and  the  assumption 
of  taking  average  values  for  gas  density  and  for  transport  coeficicnts  regardless  of 
mixture  composition,  which  is  normally  admited  in  droplet  combustion  studies,  may 
introduce  important  errors  for  the  case  of  hydrogen  and  it  is  shown  in  the  paper. 
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